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Transport properties, i.e., conductivities cr M „, Hall constant Rh, and thermopower S are studied 
for the flux state with the gauge flux cj> per plaquett, which may model the underdoped cuprates, with 
the emphasis on the particle-hole and parity/chiral symmetries. This model is reduced to the Dirac 
fermions in (2+l)D with a mass gap introduced by the antiferromagnetic (AF) long range order 
and/or the stripe formation. Without the mass gap, the Hall constant Rh and the thermopower 

S obey the two-parameter scaling laws, Rh = j^fn H (j^f . TF^r) and S — jS^fs (l^T; Tfc^r) > 
with a being the lattice constant, x the hole concentration, r the transport lifetime. The Rh and 
S show the strong temperature dependence due to the recovery of the particle-hole symmetry at 
high temperature. The ^-dependences of a xx {cc y/x) and a xy (independent of x) are in a sharp 
contradiction with the experiments. Therefore there is no signature of the particle-hole symmetry 
or the massless Dirac fermions in the underdoped cuprates even above the Neel temperature Tn- 
With the mass gap introduced by the AF order, there occurs the parity anomaly for each of the 
Dirac fermions. However the contributions from different valleys and spins cancel with each other 
to result in no spontaneous Hall effect even if the time-reversal symmetry is broken with <j> 7^ ir. 
The effects of the stripes are also studied. The diagonal and vertical (horizontal) stripes have quite 
different influence on the transport properties. The suppression of Rh occurs at low temperature 
only when (i) both the AF order and the vertical (horizontal) stripe coexist, and (ii) the average over 
the in-plane direction is taken. Discussions on the recent experiments are given from the viewpoint 
of these theoretical results. 

PACS numbers: 74.25.Fy, 74.72.-h, 71.10.Fd, 72.15.Eb 



I. INTRODUCTION 



Since the discovery of high-Tc cuprates, intensive stud- 
ies have been done on the two-dimensional antiferro- 
magnets. It is now established that the ground state 
of the 2D Hciscnberg antiferromagnet with the near- 
est neighbor interaction on the square lattice shows 
the antiferromagnetic (AF) long range ordering at zero 
temperature Q], and the low energy spin excitation can 
be described in terms of the spin wave theory. How- 
ever this does not mean that the electronic state in the 
antiferromagnets is fully understood. Compared with 
the triplet channel of the two-particle correlation func- 
tions, the single particle properties such as the angle 
resolved photo-emission spectra (ARPES)pl, and the 
singlet channel correlation functions such as the charge 
transport properties [|[ f|, ||, ^, [?|, || still remain controver- 
sial. In fact, there are two different pictures for it. One is 
the conventional spin density wave (SDW) picture || with 
the wavenumber Q — (n, ir) where both the weak and 
strong coupling regions can be smoothly connected. The 
other picture is the 7r-flux state pfj| originated from the 
resonating valence bond (RVB) idea[pT|. At half filling, 
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the 7r-fiux state is equivalent to the d-wave singlet paring 
state|||, H due to the particle-hole SU(2) symmetry^. 
The Gutzwiller projected wavefunction with the d-wave 
paring and AF orders gives a better energy compared 
with that only with the AF order Also the higher 
energy continuum of the neutron scattering spectra and 
Raman scattering spectra has been analyzed in terms of 
this 7r-flux state fia] . 

The dispersion of a single hole put into this an- 
tiferromagnet is an important issue studied inten- 
sively in terms of self-consistent Born approximation Jl6[ , 
exact diagonalizationp^7| , spinon-holon bound state 
picture [Tc|, and variational method |l9|. With the t-J 
model, all the analyses give the maximum of the hole 
dispersion at k = (±7r/2, ±7t/2). This dispersion can 
be understood in terms of the 7r-fiux picture |IM ^9[ , 
which introduces the nodal Fermi points of the spinons 
at k = (±7r/2, ±7t/2) with the dispersion similar to the d- 
wave superconductors. This fits the ARPES experiments 
in the undoped cuprates 0. 

At finite doping, the slave-boson mean-field theory of 
the t-J model predicts the state with both the singlet 
RVB and AF orders for small xgfj]. The SU(2) sym- 
metry has been employed to represent the constraint 
and the underdoped pseudogap region is characterized as 
the staggered flux state with spin-charge separation pl|, 
which can be regarded as the fluctuating state between 
the d-wave pairing state and the current order state. On 
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the other hand, the staggered flux state with the elec- 
tron coordinates with the real current ordering and peri- 
odicity doubling has been proposed for the underdoped 
cuprates |22|. 

Recently there appeared several experiments on the 
Hall coefficient Rh and the thermopower S in the heavily 
underdoped cuprates, which raised the issue of particle- 
hole symmetry || 0, gj. The Rh as well as the S show 
the strong suppression below the Neel temperature Tjv 
in some YBCO samples || while it is not the case for 
other samples 0. This strongly suggests that the trans- 
port properties below Tjv is sensitive to the oxygen chain 
ordering and/or the self-organization of electrons such as 
the stripes, which depends on the annealing procedure in 
the sample preparation. Another interesting clue is that 
in LSCO the suppression of the Hall effect is observed 
only in the vertical (horizontal) stripe state while is not 
in the diagonal stripe state for x < 0.05 ]23|, [24] . Above 
Tjv, Rh is a decreasing function of temperature, which 
remains one of the most puzzling features in the normal 
state properties H |l ^| ^q]. 

In the SDW picture, the large metallic Fermi surface 
enclosing the area of 1 — x (x: hole concentration) is 
recovered above Tjy. Therefore we expect the drastic 
change of the Hall constant from Rh = — a 2 /(l — x) to 
Rh = a 2 /x when the temperature T is lowered across 
TV- (Here a is the lattice constant.) Correspondingly 
the resistivity should be affected by the onset of the AF 
order. This is in sharp contradiction with the experi- 
ments. Even above T/v, the system behaves as the doped 
Mott insulator with the small number of holes |3[ |5[ |), 
and even a slight change of the resistivity is not observed 
at T N §. 

Considering all these clues, it is worthwhile to study 
the transport properties of the slightly doped 7r-flux state 
both above and below Tjv, which we undertake in this 
paper. We will neglect the interactions between elec- 
trons and the disorder potentials. The former is justified 
because the quasi-particle number is small at low tem- 
perature due to the reduced density of states and the 
el-el interaction is irrelevant. The latter becomes impor- 
tant at low temperature for small x where the resistivity 
shows an upturn, but is irrelevant for the temperature 
and x-range of our interest. We also study the effect of 
the stripe formation on the Hall constant and the ther- 
mopower, and show that oxygen chain ordering is crucial 
for these quantities. The deviation of the flux from tt, 
which breaks the time-reversal symmetry and produces 
current ordering, is also studied. This issue is closely re- 
lated to the parity anomaly |fj6| in (2+l)D because the 
two species of the Dirac fermions acquires the mass gap 
due to the AF ordering. Therefore the undoped and 
underdoped cuprates offers an interesting laboratory to 
study the transport properties of the Dirac fermions with 
nontrivial topological nature influenced by the AF order 
and/or the stripe formation. The Dirac fermion has been 



studied also in the context of the nodal quasi-particle in 
the d-wave supercouductors|27 . In this case the Fermi 



energy is always at E = and the particles-hole sym- 
metry remains. Under the external magnetic field H, 
the formation of the vortex lattice is crucial, which in- 
troduces the magnetic length scale Ih ~ H^ 1 ' 2 . In the 
present case, on the other hand, there occurs no Meissner 
effect, and one can see the response of the Dirac fermions 
in the uniform state, while the particle-hole symmetry is 
broken due to the shift of the Fermi energy from E = 0. 
This introduces the length scale £ x ~ x~ x l 2 , i.e., the 
inter-hole distance, and we can have the new kind of 
scaling when neither AF order nor the stripe formation 
occur. 

The plan of this paper follows. In sec. II, our model 
Hamiltonian is introduced and its spectrum under the 
external magnetic field is reviewed. In sec. [II, the elec- 
tromagnetic and thermal responses are studied for the 
Dirac fermions without stripes. The effects of the stripes 
are studied in sec. [V , and sec. [v] is devoted to discussions 
in comparison with experiments. 



II. HAMILTONIAN 

We start with the most generic Hamiltonian in the 
staggered flux state with anti-ferromagnetic and stripe 
(quasi-)order which are treated as on-site potentials in 
the mean field theory. Hereafter, we take the units in 
which h = c = 1. 

(r,r'), a 

- ^2 \} 1 + Ucr C0S (Q -r) + v cos(Q v ■ f) ci a cp a (l) 



The transfer integral if±£,r = tp±y ? ~ ^ e ' 4 (with 
r x + r y = even) represents the staggered flux <f> for each 
plaquette. In the RVB picture, cj) is generated by the su- 
perexchange interaction J and equals tt in the undoped 
case. Therefore, in our representation, the transfer inte- 
gral t is estimated as t = 0.8J, and should not be con- 
fused with t in the t-J model |I8|. In the second term, 
u a is the mean field potential for the a spin electrons 
with the wavenumber Q — (7r, tt), which originates from 
the AF order uiaf and/or the diagonal stripe formation 



^diag 



am AF + m di 



lag. 



while v is that from the 



vertical stripe formation with Q v = (tt, 0). This is the 
result of the higher order process in the stripe potential 
with the wavenumber Qstripc = (tt/M,tt/M) (diagonal 
case) or Qstripe = (Tr/M,0) (vertical case). The chemical 
potential /i is zero at the half-filling. The eigenvalues for 
each k are 
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e*(k) = ± 



4t 2 (cos 2 k x + cos 2 k y ) 



v 2 ±2\ 



At 2 cos — cos k x cos k y 



r 



v 2 {At 2 cos 2 k y + u 2 ) 



(2) 



First, we consider the case without the vertical stripe. 
In this case, the low energy electronic states are described 
in terms of the two Dirac fermions at k\ = (jr/2, tt/2) and 
&2 = ( 7r /2, — tt/2), which we call 1 and 2 sector, respec- 
tively. Measuring the momenta from these fci^-points, 
and in the continuum approximation the Hamiltonian is 
given as 



HcS = T,f d2r *lWHi*1>i*(f), 



(3) 



where ip^J (r) is the annihilation (creation) operator of 
i-sector (i — 1,2) with tr-spin, 



u a - \i 2ta(p x + e l ip y ) 
2ta(p x + e l zp y ) —u a - [i 



(4) 



and H 2 a = Ti-ic 



Here p M = — id^ and a is the 



lattice constant. As is evident from the above Hamilto- 
nian, each Dirac fermion has the mass term with posi- 
tive or negative sign, and shows the parity anomaly [p6| . 
When the external electromagnetic (EM) field is 
coupled to each Dirac fermion, the Chern-Simons term 
e^ v * A^d u A\ is generated. Therefore the Dirac fermion 
with positive (negative) u a has the right (left) chiral- 
ity. Because there are two right (R) and left (L) Dirac 
fermions, there occurs the cancellation of the Chern- 
Simons terms, and no spontaneous Hall effect results. 
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This remains true even when the flux <j> is different from 
7r and the time-reversal symmetry is broken to produce 
the current order, as long as diagonal stripes are absent, 
i.e. mdiag = 0. It is because the current pattern is stag- 
gered, and does not affect the uniform response in an 
essential way. 

Now we briefly review on the Dirac fermions in the 
presence of the external uniform magnetic field B, which 
can be analytically solved |2f| p9| . In the effective theory 
represented by eq. ([}]), the first-quantized Hamiltonian 
eq. (||) is rewritten by the replacement, p^ — > fr^ = p^ — 

eAfj,(r), where [V x A(r)] z = B. Then we define the 
following bosonic operators as 



1 



i 



~ +e zSgn(eB )$~ 



(5) 



% x + e -isgn( eB ,)|~ 1 _ (6) 



These operators satisfy the commutation relation 
[a0,ajj = 1. Using these operators and the Pauli ma- 
trices n, (i = 1, 2, 3), eq. dJ) is rewritten as 



K. 



< B, 



(a«p + aj,) Ti - isgn{eB^){a ! f > - a J,) r 2 



-U a T 3 - fll. 



(7) 



The eigenvector (n > 1) is 



\na±) = 



e n ± sgn(eB \n) ® |f) ± J e n =F sgi^e-B^iv \n - 1) <g> ||) 

I 



(8) 



with the eigenvalue ±e ncr being given by 



Kr 



St 2 a 2 \eB 6 



Ba 



B sin ■ 



(9) 
(10) 



Here \n) is the eigenvector of a^a^ with the eigenvalue 

n, |f(I)) is the eigenvector of sgn(e£?0)T 3 with the eigen- 
value +1(— 1), and the quantum index for the intra-level 
orbitals is omitted. |t(4)) should not be confused with 
a real-spin eigenvector \a), but it comes from the two- 



component nature of a Dirac fermion. Especially, the 
zero mode is given by 



|0a±) 



1 



[l±Bgn(« <r eB^)]|0)®|T)®|<T). (11) 



Therefore, there exists only one zero-mode |0<7+) or 
|0(T— ) for each sector and its energy is eo a = \u a \ or 
— | Wo- 1 depending on the direction of the external mag- 
netic field B and the chirality of Dirac fermions. We 
can also define effective current-operators as follows, 
j x = 2ta T\ and j y — 2ta(cos ^ n + sin ^ r 2 ) 
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III. ELECTROMAGNETIC AND THERMAL 
RESPONSES 

The electromagnetic and thermal linear response func- 
tions are obtained by the Kubo formula. We put the 
above solutions into the Kubo formula for a^(io) and 
approximate the effect of the relaxation by replacing lo 
by lu + i/r with the lifetime r. This approximation re- 
produces the Drudc formula for fj^ v {uj) in the simplest 
case. The contribution from 1-sector with cr-spin is pre- 
sented in the appendix |a|. The total conductivity is given 
by summing up the contribution from all sectors. This 
procedure makes anomalous terms, which come from the 
inter-band effect, cancel out with each other, and lead 
to the consistent result in the AF ordered state, i.e. no 
parity sym metr y brea king in the limit B — > 0. 

The eqs. (A5) and ( |A6| ) are valid for general cases in- 
cluding both the quantum limit (t — > oo with hnite B) 
and the semi-classica l lim it ( K eB t <C 1 <C |/j.|-r) . In 
the former case, eqs. ( [A5| ) and ( |A6| ) represent the integer 
quantum Hall effect, while in the latter case they cor- 
respond to the usual Hall effect. In the latter case, we 
obtain the same formula with that given by the Boltz- 
mann equation. For example, with only the AF order 
and at low temperature, eqs. (A5) and (lAff) lead to 



tained in a similar way as 



2^ 



8irx t 2 r 



Attx t 2 sin ^ 

(8TTXt 2 T) 2 



(12) 



l AF 



>.iy 



2?r Anxt 2 sin 



l AF 



a2 \ e \ B ■ 2 4> no , 

— sin -(13 

2irx 2 



and the Hall coefficient R B is given by 

„2 



Rh — 



a 



—— sin — . 

\e\x 2 



(14) 



independent of the presence/ absence of the AF order. 
This is reduced to the conventional result Rh = a 2 /(|e|a;) 
for the hole system with concentration x when we put 
<f> = 7r. However the x-dependences of a xx and a xy are pe- 
culiar in the absence of the AF order, namely a xx oc ^fx 
and a xy being independent of x when we neglect the IE- 
dependence of r. These are understood as follows. The 
<i xx is proportional to the density of states at the Fermi 
energy, which is proportional to *Jx for the massless Dirac 
fcrmion. As for a xy , on the other hand, the geometric 
interpretation is useful[^0|. The Hall conductivity is de- 
termined by the scattering path length, l(k) — rVe(fc), 
which is independent of the magnitude |fc| for linear dis- 
persion and hence a xy is independent of x. In eq. (A£), 



this behavior is the consequence of the dominant con- 
tribution from the zero mode of the Dirac fermions to 
a xy . On the other hand, when \f^nxt <C |7tiaf|, the 
usual dispersion e(fc) cx k 2 is relevant and the usual x- 
dependences, a xx oc x and a xy oc x result. The ther- 
mopower S at sufficiently low temperature is also ob- 



c _ k B tt 2 
e 3 



4irx f 



2m\ F 



4irx t 2 



sin ■ 



4:TTX t 2 



- m 



AF 



(15) 



The x-dependence of S is also peculiar in the absence of 
the AF order, i.e., S oc 1/^/x. 

Next we consider the temperature dependences of Rh 
and 5, which originates from the Fermi distribution 
function /^(T, e, /i(T)). Here the energy is averaged 
over ~ fcsT, and both the particle and hole branches 
contribute with the opposite signs to Rh and S when 
ksT > \fi(T)\. Furthermore, |/i(T)| is a decreasingfunc- 
tion of temperature as presented in the appendix |B|. For 
example, when vtiaf — 0, it behaves at k B T 0(t^/x) 



/ 4ttx t 2 


. 
sin — 




2 



and at k B T > 0(ty/x) 



irx t 
In 2 



-(k B T) 2 , (16) 



(17) 



k B T 



Therefore, as long as the flux order is present, the 
particle-hole symmetry is approximately restored at high 
temperature even in the doped system, and both Rh and 
S are reduced as shown in Figs. |l|(a) and |l|(b), which 
are the results of the original lattice-model. The peak 
in Rh occurs because the suppression of a xx as increas- 
ing T dominates at low temperature. After all, the low 
carrier density and the narrow gap result in these fea- 
tures, i.e., having a peak at moderate temperature and 
decreasing at high temperature. It is worthwhile to note 
that the temperature dependence of S resembles those of 
rare-earth compounds, such as YbCuAl, which have low 
carrier densities and narrow gaps [3l| . 

In the case of ttiaf = and the limit B — > 0, the fol- 
lowing scaling laws are expected by the continuum model 



Rh 
S 



a 2 [ U/x 1 
\e\x* RH \k B T' rk B T 

k B ^ / ty/x 1 

|e| \k B T' rk B T 



(18) 
(19) 



where fn H and fs are dimensionless functions. (The de- 
tails are given in the appendix ^|.) The dependence on 
t^fx/(k B T) is because of n/{k B T) being a dimensionless 
function of t^fx/(k B T), and the dependence on rk B T 
comes from the transition between the highest valence 
band and the lowest conduction band. The latter is ne- 
glected in the Boltzmann transport theory. In a physi- 
cal sense, the contribution from transitions between two 
bands would be very small when the condition 1 <C \ijl\t 
is satisfied, i.e. in the semi-classical limit. In this case, 
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FIG. 1: Temperature dependences of (a) the Hall coefficient Rh in the unit a 2 /\e\, and (b) the thermopower S in the unit 
fcs/|e| for niAF = v = 0. 




FIG. 2: Scaling of temperature dependences of (a) the Hall coefficient Rh divided by Rho = a 2 / (|e|a;), and (b) the thermopower 
S in the unit fcs/Jej for uiaf = v = 0. The open circles are scaling functions calculated in the continuum model. 



we get the one-parameter scaling lows, 

Therefore, both xRh and S are expected to scale as func- 
tions of ty/x/(kBT). Figs. ||(a) and ||(b) show the scal- 
ing behaviors of Rh and S given in Figs. |(a) and 0(b). 
We can see that the single-parameter scaling law of Rh 
works fairly well in broad temperature range. On the 
other hand, for S, it works only at low temperature. This 
is because S contains the energy integral which has one 
more energy-dimension than the integral for Rh, and S 
is more sensitive to the higher energy region where the 
lattice structure is relevant. 



In Figs. ||(a) and ||(b), the temperature dependences 
of Rh and S are shown in the presence of the AF or- 
der and the vertical (horizontal) stripe. These results 
are obtained for the lattice model in terms of the Boltz- 
mann transport theory where the inter-band effect is ne- 
glected. We can observe that the AF order (without the 
stripe formation) results in the enhancement of Rh and 
S. One reason of this enhancement is that the AF or- 
der suppresses the recovery of the particle-hole symme- 
try at high temperature. (The details are given in the 
appendix |b|.) Another reason is that the heavier mass 
rriAF makes the conductivity smaller as seen in eqs. ( p"^ ) 
and (13). The former reason is crucial both for Rh and 
S. On the other hand, the latter would be crucial only 
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for S as expected from the comparison of eqs. and 
( lilf ) . This is because the explicit dependence of a xx and 
a xy on ttiaf nearly cancel each other in Rh- Therefore, 
the enhancement of S is more drastic than that of Rh- 
As for the effect of the stripe, we will consider it closely 
in the next section. 



IV. EFFECT OF STRIPE FORMATION 

The quasi-one dimensional spin/charge ordering oc- 
curs in some cupratesjp2|, and consequently affects the 
transport properties §|) ^3). Therefore it is worth- 
while to study the effect of the stripe formation on the 
flux state. Some works have been done assuming the 
one-dimensionality, which corresponds the limit of strong 
stripe potential |33| . In this section, we give an alterna- 
tive and complemental study starting from the 2D flux 
state. It is easy to see that the effects of the stripes 
are essentially different between the diagonal and vertical 
(horizontal) ones, because the vertical (horizontal) stripe 
introduces the off-diagonal matrix elements between the 
two Dirac fermions around fa and fa, while the diagonal 
one does not. Therefore the effect of the diagonal stripe is 
similar to that of the AF order and we do not expect any 
drastic change of Rh in the case of the diagonal stripe, 
although it makes the mass gaps of different spins unbal- 
anced and leads to the modification of S. It is noted that 
the change of the mass gap modifies S as in eq. (p]) . (The 
unbalance of the mass gaps of different spins also leads to 
the ferrimagnetism.) After all, the drastic modifications 
of both R H and S are expected only when the vertical 
(horizontal) stripe formation occurs. We will focus on 
this case below. 

In the presence of the stripe, there appears the 
anisotropy in the plane. A naive expectation is that 
the ID nature of the transport along the stripe reduces 
Rh- However, this is not the case, because the stripe 
does suppress not only a xy but also one of a xx and a yy 
with Rh = <J xy l '(Ba xx a yy ) being unchanged. This is 
the case also in our explicit calculation showing that 
Rh = a /(\e\x) is a robust feature at low temperature. 
At this point, we must consider the configuration of the 
stripe in cuprates. When there are stripes, their direc- 
tion would be different in each CuC>2 layer as in LSCO, 
or there are domains of stripes with different directions 
in a layer as in twined YBCO. In the former system, 
we simply assume that the vertical and horizontal types 
occurs alternately. Then it is reasonable to average the 
contributions from different layers, because they can be 
regarded as a parallel circuit, i.e. we should take an aver- 
age of conductivity tensor not of resistivity tensor. Also 
in the latter system, the more conductive regions perco- 
late and hence dominate the conductivity of the system. 
Therefore it is reasonable to take average of a xx and a yy . 
Then the observed Rh is suppressed because only a xy is 
reduced considerably while a = (a a 
R H = cr xy /(Ba 2 ) is suppressed. 



-a yy )/2 is not, and 



Before discussing general cases, we consider some spe- 
cial cases where Rh and S can be analytically evalu- 
ated at low temperature by the Boltzmann equation. 
The first limiting condition is that IttiafI is sufficiently 
larger than the kinetic energy 0(t^/x). The second con- 



dition is \v\ <C IjtoafI or y/^t 2 



m 



AF 



< \v\. With 



these conditions, the low-energy physics would be ap- 
proximately described by quasi-particles of linear disper- 
sions at k = (ir/2, 7r/2) and (tt/2, 0), respectively. It is 
noted here that these quasi-particles are not simple two- 
component Dirac fermions, because the vertical stripe 
formation mixes the two-component Dirac fermions with 
different chiralities as mentioned above. For the case 
\v\ < \v c \ ~ 2nxt 2 jvriAF, the analysis is complicated be- 
cause we should consider two bands doped in different 
way. It is noted here that there are two upper bands and 
two lower bands for each spin degree of freedom when 
ht-af 7^ and » / 0. The analysis of Rh in this case 
is given in the appendix [e]. When |u| is larger than the 
critical value |u c |, only the second band is doped. There- 
fore, for |k c | < \v\ < \mAF\ or y/ At 2 + m 2 AF <C M, the 
result is expressed as 



Rh 



4t 2 t 2 



\e\x (t 2 +t 2 ) 2 



(21) 



where t = ty/l — \v/ttlaf\ for \v c \ < \v\ <C Im^i?!, or 

t = ty]\v\/y/W~+ m\ F -l for y/At 2 + m 2 AF < |«|. 

On the other hand, the thcrmopower ~ 
(3{Jn JfD^ 1 {Ju^a) •> where J e is the electric current den- 
sity and is the heat current density, along the 
direction fi = x,y remains almost isotropic because 
the anisotropy of the correlation functions (J^J^) and 
(JuJfH) cancels each other in the numerator and denomi- 
nator, respectively. However, S is rather sensitive to the 
change of the gap and the electronic dispersion, because 
(JuJ9) contains additional dimension of energy and is 
suppressed by the coexistence of the AF order and ver- 
tical (horizontal) stripe. For the case |v c | < \v\ <§C |?71af| 
or yj 4t 2 + m 2 AF <C \v\, it is given at low temperature by 

8ttx tt + 2m 2 



e 3 



: tt [8ttx tt + m 2 ] 



where t — iyl — \v/itiaf\ and 



\m A F\ 



(22) 



\v\ for 



|« c | < M < \m AF \, or i = t\J \v\/ y/W+m%p - 1 and 

< \v\. Com- 



At 2 



AF 



for yjAt 2 + • 



paring eq. ( |15| ) in the case<^> = 7r and eq.(22), we can see 
that the weak stripe, i.e. \v\ ~ \v c \, strongly suppresses 
the thermopower with sufficiently strong AF order, i.e., 
large |mAF|. Therefore a crucial test of this scenario is 
to measure the thermopower S in the untwined sample 
where Rh is not suppressed. 

Now let us turn to the results for Rh and S in general 
cases. Although the analytic results for Rh (eq. (pl|)) 
and S (eq. (p2|)) are valid only in the limiting cases, they 
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(a) Hall coefficient (b) Thermopower 




0.5 1 0.5 1 

k R T/t k R T/t 



FIG. 3: Temperature dependences of (a) the Hall coefficient Rh in the unit a 2 /\e\, and (b) the thermopower S in the unit 
/xV/K for {m A F,v) = (0,0), (t,0) and (t,t). 



(a) R„/R H0 at k B T = 0.02 t (b) S[k B /\e\] at k B T = 0.02 t 




FIG. 4: (a) Hall coefficient Rh divided by Rho = <a 2 / ( | e | £c) , and (b) thermopower S in the unit fes/|e|. Symbols (I), (II), and 
(III) correspond to Figs. |(I), |(II), and |(III). 



show what are crucial for each quantity. As for Rh, the 
anisotropy of the transfer integral, i.e., the anisotropy of 
the velocity, is crucial with the help of the averaging of 
the conductivity tensor. On the other hand, as for S, 
the mass reduction, i.e., the reduction of the band gap, 
is also crucial as well as the anisotropy of the transfer in- 
tegral. Therefore, interesting results are expected in the 
nontrivial case |m^f| < |u| < \/4i 2 + m 2 AF , because the 
energy dispersion is highly anisotropic and there is no 
band gap between the two middle bands. We employ the 
Boltzmann transport theory and numerically evaluated 
the correlation functions. Figs. |](a) and ||(b) shows Rh 
and S as a function of the strength of the AF order tuaf 
and the vertical stripe v, respectively. The temperature 



is fixed at Q.02t/kB- As mentioned above, the region near 
the left axis in Fig. |I](b) shows that S is once suppressed 
by the weak stripe, i.e., \v\ ~ |u c |. Then both are remark- 
ably suppressed when tuaf and v are comparable with 
each other, i.e., near the diagonal line. Figs. ||(f), ||(n) 
and ^|(III) show the equal-energy contours and Fermi sur- 
faces for the 3 points in Figs, ^(a) and ||(b). Especially, 
Fig. [H](n) belongs to the nontrivial case and shows the pe- 
culiar shape of the Fermi surface. According to ref. Q , 
a xy could be reduced in the nontrivial Fermi surface com- 
pared to that of the circular Fermi surface surrounding 
the same size of area in fc-space. Especially when the 
Fermi surface have both parts of positive convexity and 
negative one, this reduction would be very effective as 
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(I) $ = 71, m AF = 2.5t, v = 0.5t (II) § = n, m AF = 2.5t, v = 2.7t (III) § = m AF = 0.5t, v = 3.0 1 
[0, 7C] [7t, 7C] [0, 7t] [71, 7C] 




[It, 0] 





[It, 0] 




FIG. 5: Schematic view of the energy contours for three parameter sets. Thick lines in each figure represent the Fermi surface 
for x — 0.02. The reduced Brillouin zone is the square with the corners [±7r/2, ±7r/2]. Copies are also drawn over the original 
Brillouin zone for (j> — tuaf = v = 0. 



seen in the region |u| ~ |toaf|- 

Finally, the effects of the AF order and the vertical 
(horizontal) stripe on the temperature dependences of 
Rh and S are shown in Figs. ||(a) and |](b) . The AF order 
enhances the both quantities at the almost all tempera- 
ture range except for the region near the peak of Rh , and 
shift the peaks to higher temperature. This is because 
the AF order introduces the additional energy scale, i.e., 
the band gap. Then the vertical (horizontal) stripe re- 
markably suppresses both of Rh and S in the whole tem- 
perature range when v is comparable with tuaf- These 
results are parallel to the above results where the tem- 
perature is fixed. 



V. DISCUSSION 

Now we discuss the above results in comparison with 
the experiments. The most important issue is whether 
the high temperature phase above TV and the pseudo- 
gap region can be described in terms of the massless Dirac 
fermions by putting tuaf — 0. The temperature and IE- 
dependences of Rh and S appear to be qualitatively con- 
sistent with the experiments. However, when one looks at 
the x-dependence of the conductivity a xx at T — 300K 
above TV, it fits better with cx x 3 ^ 2 rather than oc x 
due to the ^-dependence of the hole mobility^. On the 



other hand, <j x 



oc 



and 



'xy 



c° in the our model, 



which is contradicting with the above experiments. Fur- 
thermore, the onset of the AF order at TV does not af- 
fect the conductivity a xx in the experiment 0, which is 
more difficult to understand from any mean field pic- 
ture. This is related to the interpretation of the insu- 
lating gap at x — 0. In the mean field picture, it is 
due to the AF ordering while there remains a large gap 
even without the AF order in the Mott insulator pic- 
ture. If the large gap disappears and crossing of the two 
bands occurs above TV, the upper band becomes relevant 



at around T x = J^fx and particle-hole symmetry will 
be recovered above T x . This massless Dirac spectrum 
also gives the novel x-dependence as described above. 
Therefore it seems that the x-dependence observed in 
experiments suggests that the Mott insulating picture is 
more appropriate for the high temperature phase. In 
other words, the Dirac fermion without the AF gap is 
never relevant to the underdoped cuprates. This is also 
consistent with the asymmetry of ARPES between hole- 
dopedQ and electron-doped [Q cuprates. Recent exper- 
iments on NCCO[Q strongly suggest that the minima of 
the electron dispersion are at k = (tt,0) and (0, 7r), and 
the particle-hole symmetry is broken. This appears to be 
consistent with the SDW picture with appropriate longer 
range hopping integrals t', and t" . However, the sign of i' 
and t" is reversed when one consider the t-t'-t"-J model 
for electron doped case, and the minimum at (tt, 0) is 
recovered by self-consistent Born approximation jl6|, and 
its results can be interpreted as the flux state with the AF 
order |}5). Therefore, it makes sense to consider the flux 
state together with the AF order, but does not without 
it. 

The AF ordered state, on the other hand, can be well 
described in terms of the mean field state with the flux 
order. In this AF state, the conventional behavior of the 
doped carriers are expected without the stripe formation. 
Therefore the effects of the stripes are the most interest- 
ing issue. It is found that Rh = a 2 /(|e|x) is a rather 
robust feature at low temperature. The suppression oc- 
curs only when the AF order and the vertical (horizontal) 
stripe coexist, and the directional average is taken within 
the plane. The thermopower is also suppressed by the 
stripe formation, and does not need the directional av- 
erage. The difference between ref. J| and ref. [[?] seems 
to be due to the sample preparation. The longer time 
annealing has been done for the former case, while the 
sample is quenched in the latter case. Therefore it is 
expected that the chain is more ordered in the former 
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case, which might help the stripe formation. The 4-fold 
symmetry observed in the magnetoresistance[[56] is inter- 
preted to be due to the induced stripe in terms of the 
in-plane external magnetic field. 

In conclusions, we have studied the transport proper- 
ties of the flux state as a model for underdoped cuprates. 
This model shows several remarkable features such as 
parity anomaly, scaling laws for Rh and S, the recovery 
of particle-hole symmetry at high temperature. Com- 
pared with the existing experimental data, in particularly 
the x-dependence of a xx , it is unlikely that this model 
describes the underdoped cuprates without the AF or- 
der. However, the flux state with the AF order, which 
gives the mass gap to the Dirac fermions, describes well 
the ordered state. In this case, the stripe order affects 
the transport properties in a nontrivial way, and we have 
discussed the experiments from this view point. 
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APPENDIX A: CONDUCTIVITY IN THE 
CONTINUUM MODEL 

This appendix is devoted to present the uniform con- 
ductivity (T^^Lo) of the continuum effective theory con- 



structed in sec. II. When the effect of the relaxation is 
approximated by replacing w by O = w + i/r with the 
life time r and we can diagonalize a given Hamiltonian, 
the linear response theory generally gives the following 
representation for the uniform conductivity. 



a,P 



a\J u \l3)(P\J^\a) 



ft + ep - e c 



(Al) 



where a and (3 are the quantum indices of eigen states, 
J M represents /z-component of the current density, and 
— fx) is the Fermi distribution function. We apply 
the above formula to the effective theory in sec. II, i.e. 
(2+l)D Dirac fermions in the external magnetic field. 
For example, in the case of the contribution from 1-sector 
with (7-spin, the following replacements are sufficient. 



E 

a,j3 
Jx ~ 



\eB\ ~ 

9tt 



E > 

n,n'— s,s'=i 



\na±) 



j x = 2ta t\ , 



j y = 2ta(cos -ti + sin -r 2 ). 



(A2) 

Here |n<r±) and e„ are defined by eqs. (|sT)-(|ll"l) in sec. II. 
Ti are the Pauli matrices in the space of |t(4))- [See 
eq. (^).] Then the contribution from 1-sector with cr-spin 
is represented by 
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i e 



\eB\ 
2tt 



E E 

n,n'=0 s,s f — ± 



/^(s'env - m) - fF(se na - /z) (n<rs|j I/ |n'c7s')(n / c7s'|j, t |na-s) 



leB[ 

2tt ^ ^ 



1 (ncrslj^n'er - s){n' a - s\j a \nas) 



^n'a i ^na 



+ 



n,n'— s=± 

fF(tn'<r ~ s/j,) - fF{e na - sfjb) (nas\j v \n' as){n' as\j^\nas) 

If {^n' a + sfj.) - /F(e„ CT - s/j.) (nas\j v \n'a - s)(n' a -s\j n \n<rs) 



i e 



n,n'=0 s=± 



^ 2 - (e n v + e«<r) 2 



(— ^(nCTslj^ln'o- — s){n' a — slj^na. 



(n<7s| i yv|n'<7 — s)(n' a — s\j a \nas) 



ZfFjfna - SfJ.) 
" y^n'a c na J 

f (— zO)3fJ (n<7s|j ly |n / (7s)(n / <7s|j M |n(7s) 



(n a s | j„ | n' a s) (nf a s \ j n \ n a s) 



(A3) 



It is noted that all distribution functions in the first line 
represent those of electrons. In the second line, we have 
done the transformation, /f(— e — m) — * 1 — /f(£ + m)> f° r 
negative energy modes in order to change the distribution 
function of the valence band to that of the quasi-hole. In 
the third line, the difference of distribution functions is 



pulled apart to each distribution function. This proce- 
dure is justified because (n a s\j ^|n' a s') ^ only for 
n' = n ± 1 as we will see below and the distribution func- 
tion ^Fi^na ±m) decreases sufficiently fast as n increases. 

The final procedure is to substitute the matrix ele- 
ments by the following explicit forms, 



(na ± |ri|n'<7±) 

(na ± |ti \n' a =p) 

(ncr ± |r2|n' cr±) 

(ner ± |r 2 |n'cr=F) 



= ±- 



T 2 



' E na -±E n+ i crT 

ina €n+la 

' E na ±E n+ i a ± 

ina^n+la 



On+l,n' + \ <Vn' + l 



+ 1,"' 



' £'n / + lcr=F-5'n / CT^ 



+ 1 



T«sgn( eJ B )- 



±isgn(e J B ) : 



ina^n+la 

I E na ±E n+ icr± 
ina^n+la 



S n +l,n' 



in' + lain' a 

E n > + i aT E n > aT 

in' + la^n'a 



&n,n' 



+1 



+1 



(A4) 



where E na ± = e na ± sgn(e£?0)u CT . (The absence of one for the dc conductivity, the contribution from 1-sector 
of zero modes is appropriately represented by E na ±.) As with cr-spin is given by 
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s=± 



, i v |u<t|[1 + sgn(su (T e J B0)] 

fF{\u.\-Bri +T - 2)2+4T - 2u 2 



t- 2 (t- 2 + 4e 2 )(e„ + ^e" 1 ) + K 4 eB (3e n - u 2 ^ 1 ) - AK 2 eB T- 2 \u a \sgn( S u a eB^ 



p 2 -1 K 2 °° 



'2n 2 ^ (e„ + £„+i)[t 2 + (e„ + e„+i) 5 



- sgn( e B) ^ 



s , . ,(^B 4 +r- 2 )[l + sgn( SUCT eB )] 

o MM - S MJ- 



r" 4 ) 2 + 8(^ + r-)r-^ + (4r"VJ 



-2^2 



4\ -22 



2^(2 



-2,2 \2 



n=l 



X 



,2 

2tt 



KbAKb, - + 4r- 2 ul) + r- 2 (Kt B ^ - t- 4 + 4r- 2 e 2 )| Ua | e ,- 1 sgn( SU(Te i? ) 

w CT (e n + £n+i) 



sgn [ sin — ] — ^— - 2^ 



2/ 2 ^e„e n+ i[T 2 + (e„ + e„+i) 2 



(A5) 



(A6) 



The terms without the distribution function represent 
the parts of the inter-band effect which remain even in 
the case \i = and UbT = 0. As for the con- 

tribution from the inter-band effect is negligible in the 
semi-classical limit (K eB<fi -C 1 <C |/u|t) compared to the 
contribution from the Fermi level. In this case, the Boltz- 
mann theory is good approximation. (The inter-band 
effect is neglected in the Boltzmann theory.) On the 
other hand, as for cr X y, the last two terms, which include 
the inter-band effect remaining in the case n = and 
k B T = 0, cancel out after summing up the contribution 
from all sectors. In other words, when the parity symme- 
try is breaking, i.e. the numbers of right and left Dirac 
fermions are unbalanced, we can not neglect the remain- 
ing inter-band effect. Especially in the case r = oo and 
B — > 0, the last line but one in cr 1 ,^ gives a contribution 
±1/2 in the unit e 2 /h as long as u a ^ 0, where the sign 
depends on u a and </>, and h is Planck's constant. How- 
ever, it is noted that the total a xy should take a integer 
value in the unit e 2 /h when it is quantized. Therefore, in 
the case where this inter-band effect is crucial, we must 
seriously consider the contribution from the bottom of 
valence bands and the top of conduction bonds, which 
are out of range of the continuum model. 



APPENDIX B: TEMPERATURE DEPENDENCE 
OF THE CHEMICAL POTENTIAL 



Here we consider the temperature dependence of the 
chemical potential /i in the limit B — > 0. When there is 
no stripe formation, the doping parameter of the contin- 
uum model is given by 



x = d * 

c=±s=± 

d 



d 2 p 



sfF(e c (p) - s/i) 



87rf 2 



(2^) 2 

p OO 

—tJ2 / deesf F (e-sfx)(Bl) 

in t\ „__(- J\m A F\ 



Here d s is the number of spin degrees of freedom, i.e. 
d s = 2, df is the total number of inner degrees of freedom, 
i.e. df = (left + right) x d s — 4, and 



e±(p) = -\/(2to) 2 [p 2 ± 2 cos - p x p y 



AF ' 



(B2) 



The sign of x is taken as it is positive when [i < 0. 
In the low-temperature approximation, k B T -C (|u[ 



I^afIi we can estimate the right hand side of eq. (Bl) 
by using the sharpness of the Fermi distribution function 
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as follows, 



(If 

X = IM 2 s g n (v)6(\v\ - \m A F\) 



A* 2 - m\ F + — (k B T) 2 



(B3) 



where df — 4 is substituted. It is noted that, fixing 
the doping parameter and the temperature, the above 
formula is a increasing function of \m A F |- Therefore, the 
AF order tuaf suppresses the recovery of the particle- 
hole symmetry. 



Then, in the hole doping case, i.e. /i < — IjtoafI, the 
temperature dependence of \i is given by 



/ 4:7TXt 2 



. 

sin — 

2 



m\ F ~\{k B T) 2 . (134) 



where we have used df =4. 

On the other hand, in the high-temperature approxi- 
mation , fcg T ^> \/i\, we can estimate the right hand side 
of eq. (Bl) as follows 



87rf 2 



y- 



s=± @ J f3\m A F\ 



dy 



s y 



1 -f e V-sf3ii ' 



Ant 2 



sin ■ 



P\m AF \ 

1 _|_ e /3\m AF \ 



+ ln(l + e-' 9 l m - 4F l) + --- 



(B5) 



where (3 = l/(k B T) and we have expanded the integrand 
in /3/i. Finally, in the hole doping case, i.e. /i < 0, the 
temperature dependence of p is given by 



irxt 2 f3 



sm ■ 



r g^L r + ln(l + e-/'l™^|)' 



(B6) 



APPENDIX C: SCALING LAWS OF R H AND S 



In this appendix, we present the derivation of the scal- 
ing laws of Rh and S in the case of m A p = "^diag = v = 
and in the limit B — ► 0. In order to see the tem- 
perature dependence of the Fermi distribution function 
we introduce the function f(j3£) — f F (£) where 
= l/(k B T). Then, from eq. pD) with m AF = 0, the 



doping parameter x is given by the following equation, 



2nt 2 



E 

\ " 



2nt 2 



sin 1 1 ,4± J" 



de es f(/3(e - sfj,)) 



dy ys f(y - sf3/j), 



(CI) 

It is easy to see that t^fxj (k B T) is represented by a func- 
tion of n/(k B T). Therefore, when we can consider the 
inverse of the function, p/{k B T) is represented by a func- 
tion oityfx/{k B T). 

In the same way, by making the energy integral dimcn- 
sionless, the response functions are represented as 



>xy 



2tt 

i- 

2n 



sm • 



sm — 
2 



df f 1 

V ~d~y yV ~ Sl3 ^ + ( 2 ~ ^ V ~ S/3 ^ 



2a 2 eB ,„ /2r\ 2 v-- f°° , 

— —(^) 2 (j) E L d y ; 



0*2* = 



2n 



sm ■ 









2 d y 

=± Jo 



df 



dy 



1 



(*»)J 



sy(y - s0fi) — (y - s0fi) - [ - - f(y-s(3fi) 



fa 



(*') J 



(C2) 
(C3) 
(C4) 



where a® is the response function given by the correla- cr-spin is given by 
tion function of the electric current density eJ and the , 

heat current density J®. In the first-quantized repre- = (2to) 2 \p x +cos—p y I — /j,j x , 

sentation, the heat current density for the 1-sector with V 2 

jy = (2to) 2 I Py +COS-p x \ - (ijy 



(C5) 
(C6) 
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This is the continuum version of the heat current den- 
sity derived in the lattice model, which is presented in 
the appendix [d| The right hand side of each response 
function is a function of tyfxj '(fcsT) and rfc^T, be- 
cause fi/(ksT) is a function of ty/x/(ksT). Therefore, 
xRh = x cr xy l (Ba^. x ) and S = &x X / (Ta xx ) is functions 
of tyfx/ '(ksT) and rfc^T. In the semi-classical limit, the 
dependence of rfc^T is negligible, and xRh and S scale 
as functions of t^fx/ (kBT). 



APPENDIX D: HEAT CURRENT DENSITY IN 
THE LATTICE MODEL 

The definition of the heat current density is not 
straightforward as against that of the electric current 
density. Here we present the heat current density of the 
system represented by eq. @) . The Hamiltonian is rewrit- 
ten as H = J^phf where 



h ,-• 



Vp„ 



r-\-5,(T 



-^2{V Pa - \i)c\ a cp a , 

a 

\u a cos(Q ■ + v cos(Q v • f) 



(Dl) 
(D2) 



and the sum of the unit vector 5 runs ±x and ±y. It 
is clear that hp is interpreted as a local heat density. 
The heat current density would be conceptually defined 
by the averaging the product of the local velocity and 
hp. However, it is difficult to define the local velocity, 
especially in the second quantized formalism. Here we 
define the heat current density from the analogy with 
the conceptual definition. First, using hp, we introduce 
the quantity R as follows, R — V^ 1 ^pr hp, where V is 
the volume (area) of the system. This quantity has the 
dimension of heat density multiplied by length. Then, we 
define the heat current density as a time derivative of R, 
i.e. J® — i[H, R]. From the dimension of R, J® correctly 
has the dimension of heat current density. Finally, the 
explicit form of is given by 



J h - fiJ, 



V 



£ 



<5.<5' 



(D3) 



- c f - - i 

+S,r P+S+S' ,<r 



- 6 (V-.,f + Vp) tx,i - ct - Cpa 
9 L ' K r+6 ~ > r+S,r p+S,a 



If £ £ (~4i 2 ) jl sin fc M I cos k^ + cos - cos(e^, v k v ) j ct (T c fe 

~ ■ , f <t> t 

2vt y sin fc„ < — cos — cl cp , ^ 

y v 1 4 k<J fe+Qx: 



i sm — c- Ct , a 

4 kcr k+Q y ,a 



f = 7 £ E s w i+i* c ™ = v E E 



2i /t sin fcp 



r <5,<T 



m x (pi 
cos— cl c,- + ( — lrisin— cl c,-,n 



4 fee k+Q,a 



(D4) 
(D5) 



where Q = (it, it), Q x = (tt, 0), and Q y = (0, tt). The 
symbol (-l)^ 1 means that (~l) x = 1 and (-1)» = -1. 

In a physical sense, J E is the energy current density, and 
J is the current density. 

As for the continuum limit, the following transforma- 
tion of bases are performed before taking the limit. 



(D6) 



C k+Q y ,<7 



C k+Q,cr 




~u+ o ' 






C k + Qz; ,<J 




U- 




. ^k,2o . 



where ipt . (i — 1,2) is a two-component operator re- 
spectively, and 



(D7) 



U± = ~~^( T 3 i T i) ( cos o T i± sin — r 2 



When there is no stripe, ipz ■ with k ~ (7r/2,7r/2) cor- 
responds to the Dirac fermion of i-sector (i = 1,2) with 
cr-spin. 
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APPENDIX E: R H IN THE CASE \v\ < \v c \ 

Here we present the analysis of Rh when two lower 
bands are doped. Around the tops of these two bands, 
their dispersions are approximated as 



where x+ and X- are doping parameters for the first and 
second bands respectively. Then, from the conservation 
of the total doping x and the commonness of the chemical 
potential /i, x+ and X- are determined by the following 
equations, 



e±(p) = -y/(2ta)'p* 



(2t±a) 2 p 2 y 



m±, 



(El) 



where t± = tyl ± \v/rriAF\ and rh± = \uiaf | ± \v\. In 
this case, by considering the two bands separately and 
then adding the contributions at the stage of conductiv- 
ity, we can get the Hall coefficient Rh — a 2 /(|e|a;) at 
low temperature, where x is the total doping parameter. 
However, when there are domains of stripes with differ- 
ent directions, we should take an average of conductivity 
tensor. Then, we obtain the following result. 



.2 4 



Rh — r~i 



(x. 



2x 



0- 



(E2) 



x+ 



(E3) 



35-4/1 ~ 



m AF 




m AF 



\vm AF \ 
2nt 2 



(E4) 



It is noted here that x± includes the contribution of both 
spin degrees of freedom. When |v| is larger than the 
critical value \v c \, i.e., 



\v\ > \v c \ = \m AF \ 



s/2-KXt 

mAF 



\ 



1 / V27rxi 
4 \ m AF 



\[2/KXt 

m AF 



(E5) 



only the second band is doped, i.e., x+ = and X- = x. 

We can also analyze the thermopower S by the similar 
way. However, its expression is more complicated than 
that of Rh and is not suggestive. Therefore, we do not 



give the explicit analysis, but pointed out that, in the re- 
gion < \v | < |u c L S is highly reduced when |v| increases 
as shown in Fig. 0(b). 
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